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ABSTRACT 


The thermal control system of a spacecraft is designed to maintain all spacecraft 
components within their specified operating temperature limits throughout all phases of a 
spacecraft’s mission. In order to verify and aid in such a design process, a thermal 
analysis of the system must be conducted. A thermal model of the spacecraft is used to 
simulate its behavior under given thermal environments and boundary conditions so that 
temperature predictions can be made. 

The focus of this thesis is to develop and analyze thermal models of PANSAT 
which describe its thermal behavior while it is in orbit and also prior to its insertion in its 
orbit (while it is still in the shuttle). The results of these analyses will serve to help in the 
thermal design and performance of PANSAT. This thesis completes the thermal model 
prerequisites of the STS 95 space shuttle mission hitchhiker program for PANSAT. The 
emphasis of this thesis is to develop a model that will allow the prediction of the 
temperatures of all the electrical components including the temperature sensitive 


electrical components of PANSAT such as batteries over a complete orbit of the satellite. 
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I. INTRODUCTION 


The objective of this thesis is to develop a model that simulates the transient and 
steady state thermal behavior of PANSAT (Petite Amateur Navy Satellite) while in 
flight/orbit under a given set of external thermal environments and boundary conditions. 
Once required model accuracy is attained, the model will be used to simulate different 
thermal environments in order to determine other design parameters, such as the amount 
and type of thermal insulation and shielding that is needed in order to maintain all 
spacecraft components within their specified operating temperature limits. The most 
temperature sensitive components are the batteries which have narrow temperature 
operating margins that must be maintained to ensure proper operation. 

This thesis covers two distinct phases of the mission. The first phase is the time 
interval while PANSAT is in the Shuttle (as the Shuttle orbits the Earth), awaiting 
insertion into its orbit. The second period is the in-orbit phase of the mission when 
PANSAT has been inserted into its orbit and becomes fully operational. Each of the 
aforementioned periods will be modeled separately and analyzed under simulated 
thermal environments and boundary conditions. These models will be used for steady 
state analyses and transient analyses to predict real time thermal performance of the 


spacecraft systems and individual components. 


A. PANSAT BACKGROUND 


PANSAT is a small, direct sequence spread spectrum communications satellite 
being developed by the students, staff, and faculty at the Naval Postgraduate School 
(NPS). PANSAT serves as a hands-on educational project for officer students in 
numerous academic disciplines to gain knowledge in designing, building, testing, and 
Operating small satellites. Additionally, PANSAT will serve as a proof of concept of low 
cost, spread spectrum, digital over-the horizon communications. Amateur radio operators 
will use PANSAT to store and forward communications once it is in orbit and ready for 
full operations. The Naval Postgraduate School mission control center will serve as the 
ground control center for the in-orbit phase of the mission. 

PANSAT consists of a simple design in that it has no attitude control or 
propulsion system and weighs approximately 150 pounds. Hence, PANSAT is to be a 
tumbling spacecraft. PANSAT is currently manifested on the space shuttle mission (STS 
95) which is scheduled to launch from Kennedy Space Center in Cape Canaveral, Florida 
in October 1998. PANSAT will be inserted into orbit by means of the space shuttle 
hitchhiker program ejection system. Silicon solar cells (manufactured by Spectralab Inc.) 
and nickel-cadmium battery cells (manufactured by Sanyo) will provide the power for 
PANSAT while it is in orbit. PANSAT will be in a low earth orbit with an inclination of 
approximately 28.5 degrees and an altitude of approximately 310 nautical miles (574 
km). PANSAT has a two year design lifetime which is primarily a function the 
aerodynamic drag and radiation experienced by satellites at the low earth orbit altitudes 
[Ref. 3]. The aerodynamic drag, a non-conservative force, continuously extracts energy 
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from the orbit; thereby, causing the spacecraft to spiral into the earth’s atmosphere and to 
subsequently burn up upon reentry [Ref. 3]. The high altitude of the STS-95 mission 


however, offers an orbit lifetime in excess of four years. 


B. SPACECRAFT THERMAL CONTROL 


The thermal control system of a spacecraft is designed to maintain all spacecraft 
components within their designed operating temperature limits throughout all mission 
phases. In order to complete the thermal design of a spacecraft, a thermal model of the 
spacecraft must be made which accurately depicts its thermal behavior under the given 
boundary and thermal environmental conditions. The focus of this thesis is to develop a 
thermal model that simulates the transient and steady state thermal behavior of PANSAT 
and to conduct thermal analyses to test the model and hence improve the thermal design 
and performance of PANSAT. 

This thesis covers two periods. The first period is the time interval while 
PANSAT is in the Shuttle (as the shuttle orbits the Earth), awaiting insertion into its 
orbit. The second period is the time interval after PANSAT has been inserted into its 
orbit. Each of the aforementioned periods will be modeled separately and analyzed under 
simulated thermal environments . 

1. Spacecraft Thermal Design 

The thermal design of a spacecraft is affected by the power, attitude and control, 
communications, and structural subsystems [Ref. 6]. The interaction between these 
systems necessitate tradeoffs that significantly affect the design of the thermal control 


system. Normally, small, unmanned, low earth orbit satellites can be controlled passively 


[Ref. 7]. Passive thermal control is preferred over active thermal control because of its 
simplicity, reliability, and lower costs. Active thermal contro] uses heaters, active 
mechanical, or active fluids devices in the design whereas passive thermal contro] use 
primarily thermal coatings, surface finishes, insulation, and equipment placement, and 
material selection in its design. Since PANSAT its in its advanced development stage, 
thermal coatings, surface finishes, and insulation will be considered for use at this point 
in the project. Normally, the thermal engineer interacts with the structural engineer, 
electrical power engineer, and attitude and control engineer in the early stages of the 
design in order to achieve optimum thermal design, system integration, and performance. 
The first step in the thermal design process is to determine the component 
temperature limits, internal power dissipation, and environmental heat inputs incident on 
the exterior surfaces of the spacecraft [Ref. 7]. The thermal design of a spacecraft is an 
iterative process based on the successive thermal analyses of the thermal model under 
various boundary conditions and external thermal environments. Figure | illustrates the 
various inputs into the thermal design process and the outputs of the thermal analysis of a 
given model. The results of the analysis are subsequently used to improved the design 
and verify the tentative design. This iterative process continues until the desired 


optimization and system integration is achieved. 
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Figure 1. Thermal design process (Ref. 1) 










II. THERMAL ENVIRONMENT 


A. ORBIT ANALYSIS 


As a member of the space shuttle hitchhiker program, PANSAT is considered to 
be a secondary payload; hence, it will have to accept the orbital elements of the primary 
space shuttle payload [Ref. 12]. PANSAT will be placed in a low-earth orbit of 
approximately 310 nautical miles with an inclination of approximately 28.5 degrees. It 
will be subject to various thermal environments depending on its location and attitude 
within the orbit. These inputs include ultraviolet solar radiation, ultraviolet earth reflected 
solar radiation, and infrared earth emitted radiation. The rate at which the spacecraft and 
the space environment exchange heat is primarily dependent upon the temperature of the 
two bodies. Since the spacecraft thermal environment heat varies with spacecraft altitude, 
spacecraft attitude, and spacecraft position in the orbit. The heat fluxes incident upon 
PANSAT vary also. The classical orbital elements, explained in Table 2.1 and shown in 


Figure 2, completely describe the orbit of any object. 


ORBITAL ELEMENT = |SYMBOL|  _ DESCRIPTION 


| Inclination ] Angle between earth’s equatorial plane | 
and the satellite orbital plane 


Right ascension of the ascending Q Angle from vernal equinox to 
































node 





ascending node 








Argument of perigee Angle between the ascending node and 
the point of periapsis, measured in 


satellites orbital plane 









T Angle from the point of periapsis to 


satellite’s position 


Table 2.1 Classical orbital elements 


rue anomaly 





Figure 2. Classical orbital elements 


These classical orbital elements determine the duration of the sunlight and eclipse 
portions of the orbit. The duration of the sunlight portion of the orbit and the eclipse 
portion of the orbit have a significant impact on systems operations, the electrical power 
system and consequently the thermal control system. Low earth orbit spacecraft have 
shorter and more frequent eclipse periods than geosynchronous earth orbit spacecraft 
which result in higher discharge and recharges rates for the batteries. The excess heat 
generated by the cyclic charge and discharge creates more heat for the thermal control 
system to manage. 

Figure 3 shows the period, minimum and maximum eclipse time for circular 
orbits at various inclinations. For PANSAT’S orbit (=28.5, altitude=310 nautical miles, 
e=0) the period is approximately 95 minutes with an maximum eclipse duration of 


approximately 36 minutes [Ref. 10]. 
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Figure 3. Orbital period, maximum and minimum eclipse time for circular orbits 
(Ref. 10) 


The minimum eclipse duration increases as the inclination of the spacecraft decreases for 
a given altitude. In addition, the period increases or decreases with a corresponding 


change in altitude [Ref.10]. 


B. EXTERNAL THERMAL ENVIRONMENT 


Figure 4 depicts the various heat sources and sinks that a spacecraft in low earth 
orbit experiences [Ref. 12]. The heat sources are direct solar radiation, reflected solar 
radiation from the Earth, emitted infrared radiation from the earth, and internal heat 
dissipation from the internal components of the spacecraft. The majority of the Solar and 
reflected radiation is in the ultraviolet region [Ref. 6]. The spacecraft exchanges heat 
with space, a heat sink, via radiation. The solar intensity varies approximately 3.5% due 
to the Earth’s elliptical orbit [Ref. 6]. The environmental heat inputs are defined using 
solar vector angles usually referred to as the beta angles. A beta angle of 23.5 degree is 
equivalent to a solar incidence angle of zero upon the spacecraft which is the worst case 
angle. 

The attitude of a spacecraft normally has a direct impact on most spacecraft; 
however, because PANSAT has no attitude control system, the spacecraft will tumble 
throughout its orbit. The most likely orientations will be investigated in this thesis. The 
non-tumbling orientation would have a tendency to set up thermal gradients because of 


the prolonged exposure of a particular surface to a particular environment. 
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Figure 4. Spacecraft thermal environment in LEO (Ref. 14) 
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Ill. SPACE CRAFT CONFIGURATION AND DESIGN 


A PHYSICAL DESCRIPTION OF PANSAT 


A schematic of PANSAT is shown in Figure 5. PANSAT is a 26 sided polyhedron 
and has a highly symmetric structure. This configuration was chosen because PANSAT 
has no attitude control or propulsion system. PANSAT is divided into three modular 
sections: the upper deck, the mid-deck and the lower deck. These sections are shown in 
more detail in Figures 6, 7, and 8. These sections consist of a top panel, a bottom panel, 
four upper deck and four lower deck panels, and eight mid-deck panels. The upper and 
lower sections have four triangular plates with the antennae mounted on the plates of the 
upper section. These antenna subassembly plates have testing ports and lifting eyes for 
PANSAT and have bare surfaces. The remaining panels are each covered with a layer of 
silicon solar cells manufactured by Spectralab Inc., and serve as a solar collector. These 
cells (model K6700A) have a mass of 55 mg/cm’ and an efficiency of 14.6% (based on 
the manufacturer’s data sheet). The launch vehicle interface (LVI) connects PANSAT to 
the base of the (GAS) canister via a base plate or bottom plate which is in turn connected 


to the space shuttle. 


3 


—_ 
&76 oT ee on 


= 4,78 In. 





Figure 5. PANSAT exterior structure with the launch vehicle interface 
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Figure 6. Upper deck configuration 


The configuration of the upper deck of PANSAT is shown in Figure 6. The upper 
module contains the mass storage containers A and B, the upper equipment plate, four 
triangular antenna subassembly plates, a top panel and four upper deck panels, and a 
temperature sensor multiplexer (TMUX). The mass storage compartments are assumed to 
have perfect contact with each other and the TMUX. The TMUX Is also assumed to have 
perfect contact with the upper equipment plate. The electrical power system (EPS) 1s not 
part of the upper section but is mounted on the bottom of the upper equipment plate with 
fasteners. All panels except the antenna subassembly panels are covered with silicon solar 
cells on the outside. 

The top panel is a square whose side is 7.125 in. in length and has a thickness of 
0.25 in. The upper equipment plate 1s an octagon whose sides are also 7.125 inches in 


length. Each side has a 125-135 degree taper with a thickness of 0.25 inches. The top 


BS) 


panel is attached to the other four upper deck panels via six 4-40 screws. The four upper 


deck panels are 7.125 in. x 8.250 in. Each panel has a 35 degree taper. 
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Figure 7. PANSAT mid-deck configuration 


The mid-deck section of PANSAT (shown in Figure 7) consists of eight deck 
panels, the electrical power system (EPS), battery compartments A and B, the digital 
control system (DCS), and the lower equipment plate (LEP). The radio frequency section 
is mounted on the bottom of the lower equipment plate although it is not part of the mid- 
deck section. There is a 0.5 in. gap between the DCS and the battery compartments. The 
DCS and the battery compartments are mounted on the LEP with fasteners. The exterior 
of all eight panels is covered with silicon solar cells. These panels are 7.125 in. x 7.125 
in. on a side, and 0.063 inches thick. These panels exchange heat externally by radiation 
with the walls of the GAS canister, while they primarily exchange most of their heat by 


conduction with the internal components of PANSAT . 
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Figure 8. PANSAT lower deck configuration 


The lower module of PANSAT (shown in Figure 8) contains the radio frequency 
section, the support cylinder, the base plate or bottom plate, four lower deck panels and 
four triangular plates. The support cylinder made of aluminum is hollow and attaches 
directly to the lower equipment plate (LEP) and the base plate of PANSAT. The LVI 
attaches to the base plate and connects PANSAT to the base of the GAS canister. In this 
model, the LVI is lumped together with the base of the GAS canister and together 


assumed to be thermally massive. 


Ie 





IV. THERMAL MODELING 


A. HEAT TRANSFER 
1. Conduction 


Conduction heat transfer occurs between two solid or fluid objects in physical 
contact with each other under the influence of a temperature gradient [Ref. 1]. Heat 
transfer occurs across the medium. Conduction is simply the transfer of kinetic energy 
from more energetic particles to the less energetic particles [Ref. 4] . The basic law of 


one dimensional conduction is governed by the Fourier’s equation. 


aT 
= —kA—— 4.1 
Q X [4.1] 


Where, 
Q= the rate of heat conduction along the x axis ,W 
k= thermal conductivity, a physical property of the material , W/m-K 


A= cross sectional area of the path normal to the x axis ,m” 
K 
—— = temperature gradient along the conduction path, — 
dX m 
The thermal conductivity (k) varies with temperature; but, may be assumed to be 
constant over the a given temperature range [Ref. 6]. The conduction heat transfer is 
analogous to the Ohm’s law for electricity [Ref. 1]. The heat flow ,Q, represents the 
current. The temperature gradient, dT, represents the voltage and kA/dx represents the 


conductance which is the reciprocal of the thermal resistance [Ref. 1]. The conductances, 


Lo 


the reciprocal of the thermal resistances, account for the heat transfer between nodes as 


shown in Equation 4.2 and Figure 9. 
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Figure 9. Geometry to calculate linear conductances 


C= [4.2] 


ie 

where, 

C = Linear conductance (W/K) 

k = Thermal conductivity (W/m-K) 


A = normal cross sectional area between nodes(m*) 
L =distance between nodes(center to center ) (m) 


Contact between to surfaces are rarely perfect and the pressure applied to any 
component is rarely uniformly distributed. Thus, there is a contact resistance at the 
surface interfaces of the two objects. This contact resistance plays an important role 
when boxes are mounted to surfaces within the spacecraft. (i.e., DCS, battery 
compartments, EPS, and RF COMMS). The contact conductance varies with the hardness 


of the material, surface finishes, surface contact area, and the pressure applied [Ref. 6]. 
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Gap filler material is used in order to fill the void due to surface roughness and increase 
thermal conductivity [Ref. 6]. 


2. Radiation 


Radiation is electromagnetic energy that is emitted from all bodies whose 
temperatures are above absolute zero [Ref. 7]. Radiative heat transfer is defined by the 
Stefan-Boltzmann law and varies as the fourth power of the temperature [Ref. 7]. Hence, 
at lower temperatures, larger areas are required to reject a given amount of heat. Radiative 
heat transfer does not required objects to be in physical contact. Radiation can take place 
in a vacuum and is a function of temperature of emitting and receiving bodies (at least 
One temperature must be above absolute zero) , the surface properties of the material of 
the bodies, the intervening medium, and the relative geometry. The total radiation 
incident upon a spacecraft is either absorbed, reflected, or transmitted [Ref 7]. The 
radiation that is absorbed must be managed by the thermal control system. Equation 4.3 


describes radiation heat transfer. 


q = €0T" [4.3] 
Where, 

q= total energy radiated per unit area 

€ = emissivity 

6 = Stefan Boltzmann Constant 


T= absolute Temperature in Kelvin 
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The relative geometry of the surfaces plays a major role in radiation heat transfer. 
This geometry is characterized by a view factor which relates the fraction of energy 
emitted from one surface that is intercepted by another surface [Ref. 4 ]. Figure 10 shows 


the radiation heat exchange between two surfaces and equations 4.4 and 4.5 depict the 


view factor of surface 1 with respect to surface 2 [Ref. 4]. 





Figure 10. View factor associated with the radiation heat exchange between two 
surfaces (Ref. 1) 


1 dA, cos¢, dA, cos¢ 
F,=— | {| ———-—— [4.4] 
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Where, 

F,.= view factor from surface | to surface 2 
A, = areaof surface 1 (m’) 

A» = area of surface 2 (m7) 


d= angle between normal and surface | 


d2= angle between normal and surface 2 


The amount of radiation emitted from surface 1 that is incident on surface 2 ts 


given by the following relationship. 


On = Ane Oe en) [4.5] 


where, 

Q)2= the net heat exchange between two diffuse gray surfaces land 2 (W) 
A ;= area of surface 1 (m’) 

Oo =Stefan Boltzmann constant (W/m’K*) 

T, = temperature of the surface (K) 

T2 = temperature of the surface (K) 


F,2= view factor from surface 1 to surface 2 


Ze 


3. Convection 


Convection is heat transfer between a flowing liquid and a solid interface [Ref. 
13]. Convection plays a minor role in low-Earth orbit, unmanned spacecraft and is not 


discussed in this thesis [Ref. 13]. 


B. LUMPED PARAMETER MODELING 
1. Thermal Nodes 


The thermal behavior of satellites is completely described by the thermal 
capacities (nodes), conductors, and radiative couplings which form the analytical model 
[Ref. 7]. Although the surface properties of the nodes degrade over an extended period of 
time, the optical properties (absorptivity and emissivity) change minimally for an 2-4 year 
design lifetime of PANSAT. For longer in-orbit phases, and greater altitudes, the 
absorptivity increases and the infrared emissivity remains relatively constant. Therefore, 
end of life (EOL) optical properties are used in thermal control system design. 

The size and selection of nodes are chosen so that they are isothermal. A material 
with a high thermal conductivity such as aluminum can be considered isothermal. The 


nodes contain all the mass, geometry, and optical properties. 


C. ASSUMPTIONS 


1. PANSAT In The GAS Canister 
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Figure 11. PANSAT in the gas canister 


The thermal model is based on a lumped nodal analysis of the various key 
components of this configuration [Ref. 2]. In all, the physical setup comprising PANSAT 
and the GAS canister have been broken down into 12 nodes, with each node representing 


a region which may be treated as isothermal. In the model of Figure 11, the nodes are 


ae 


allowed to exchange heat with one another by conduction depending on thermal contact, 
and/or by radiation depending on the view factor geometry [Ref.4]. A schematic of the 
nodal network is shown in Figure 12, and the details of each node are provided in 


Appendices A and B. 


The major assumptions in Setting up the model are the following: 


e The solar constant (solar radiation flux) is taken to be invariant at 1353 W/m? over the 
entire orbit [Ref. 1]. 


e The top, mid-deck and bottom sections and enclosures, each make up an isothermal 
node and have uniform radiosity and irradiation characteristics [Ref. 3]. 


e All surfaces exhibit opaque, diffuse, gray surface behavior. 


e Each node is isothermal with all the material and thermal properties of the region it 
represents being lumped into the node. 


e The mass of the base of the GAS canister coupled with the shuttle is sufficiently large 
such that its temperature does not vary over the course of this analysis. 


e Only the top section of PANSAT receives direct solar radiation in this model. The 
mid-deck and lower sections are not exposed to space while PANSAT is in the GAS 
canister due to the close proximity of their walls. 


e No single section of PANSAT sees the entire inner surface of the can - different 
sections see different portions of the GAS canister depending on their location and 
orientation. 


e The four strip heaters located in the lower part of the GAS canister see only the lower 
deck section of PANSAT, and each provides equal heat distribution to PANSAT. 


e All the electronic boxes and components inside PANSAT which are fastened or 
bolted to one another are in perfect thermal contact. 


e All of the triangular antenna subassembly panels are also taken to be covered with 
solar cells. 

e The mass of the solar cells themselves are neglected in the calculations of the nodal 
masses. 
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e The GAS canister can receive solar radiation in an unobstructed manner. 


e A representative temperature of 250 K is assumed as the initial temperature in this 
analysis. Such a choice however is not restrictive in the interpretation of the results as 
discussed later in Section 9. 


e All satellite materials are isotropic and their properties are constant over the 
temperature range of interest. 





Figure 12. The thermal network model (in the canister) 


2. PANSAT in Orbit 
PANSAT is in its orbit and fully operational in this model. The LVI is still 
attached to the spacecraft as shown in Figure 5. 


The major assumptions concerning this model are: 


Zan 


The top, mid-deck and bottom sections and enclosures, each make up an isothermal 
node and have uniform radiosity and irradiation characteristics [Ref. 3]. 


All surfaces exhibit opaque, diffuse, gray surface behavior. 


Each node is isothermal with all the material and thermal properties of the region it 
represents being lumped into the node. 


All the electronic boxes and components inside PANSAT which are fastened or 
bolted to one another are in perfect thermal contact. 


All of the triangular antenna subassembly panels are also taken to be covered with 
solar cells. 


The mass of the solar cells themselves are neglected in the calculations of the nodal 
masses. 


A representative temperature of 250 K is assumed as the initial temperature in this 
analysis. 


All satellite materials are isotropic and their properties are constant over the 
temperature range of interest. 
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Figure 13. The thermal network model (in orbit) 


D. BOUNDARY CONDITIONS 
1. PANSAT in the GAS Canister 


Nodes 2,3,4,5,9, and 10 are the interior nodes and do not exchange heat directly 
with the external ambient. Nodes 1,6,7,8,11,12 on the other hand have exterior surfaces 
that can exchange heat with the external thermal environment. Of these, only nodes 1, 11 
and 12 representing the top panels of PANSAT, and the sides of the GAS canister 
respectively, receive solar radiation. The base of the GAS canister serves effectively as a 
heat sink which tends to transfer all the heat from PANSAT into the Space Shuttle. The 
exterior surface of the GAS canister is covered with thermal beta cloth. All of the kapton 
strip heaters are taken to be either on, or off, at a particular instant in time. The battery 


compartments (nodes 9 and 10), the EPS (node 3), the RF communications box (node5), 
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and the DCS (node 4) are not powered during these storage conditions and produce no 


heat while PANSAT is in the GAS canister. 


2. PANSAT in Orbit 


The boundary conditions in this model are similar to the previous model with 
following exceptions. The GAS canister ( nodes 11 and 12) is absent. The launch 
vehicle interface (LVI) is node 11 in this model. The batteries, DCS, RF COMMS, and 
EPS are powered during this phase. The interior nodes are the same as in the previous 
phase. Nodes 1,7,8,9,6, and 11 are the exterior nodes and can exchange heat with the 
external thermal environment. 


Equation 4.2 was used to calculate the conduction heat transfer conductance 
between the various nodes in thermal contact [Ref. 1]. This value of conductance is based 
on the geometry of the section and the material used. It is not temperature dependent 
since material properties are assumed to be constant. The contact area is assumed uniform 
with an even pressure distribution. These values for the conductances are tabulated and 
listed in Appendix A. The area is the common cross sectional area of the two surfaces in 
contact. The area is the smaller of the two areas if the surfaces are of unequal sizes. The 
length, L, is the distance between the nodes normal to this cross sectional area as shown 


in Figure 10. 


Equation 4.4 refers to heat transfer by radiation and is dependent upon the current 
temperature of the nodes in question [Ref.1]. In the model, the radiation conductance 
v ‘ue is calculated for each time interval to get the most accurate value for the 
conductance at that instant in time. Thus the value of the radiation conductances are 
modified in the time integration process as the temperatures approach their steady state 


values. 
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Equations 4.3 and 4.5 are used to establish the heat balance at each node [Ref. 1]. 
The net heat exchange at each node is determined, and is set equal to zero for a steady 
state solution. The transient solution, however, 1s dependent upon the thermal capacitance 


of each node and the net energy in-flow into the node. 
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V. METHODOLOGY 


A. PROGRAM CODE 
1. PANSAT in the GAS Canister 


The models described above and the associated time integration procedure were 
implemented with the help of Matlab (the source code listing for which has been provided 
in Appendix D). The time integration procedure is carried out for each of the twelve 
individual nodes (base excluded) simultaneously in a coupled manner to obtain the 
transient values for each nodal temperature, and to study their evolution to their 
respective steady state values. 

Using the model and assumptions described in Section 4, PANSAT has been 
divided into ten different nodes (nodes 1-10) while the GAS canister was divided into 
three nodes (11,12 and the base). See Appendices A and B for the composition and 
description of each node. The twelve nodes, not including the base node, give twelve 
nodal energy balance equations at a given instant in time. These twelve equations form a 
system of 12 coupled nonlinear (due to the radiation terms) ordinary differential equations 
that are solved simultaneously using the “ode45” differential equation solver in Matlab. 
As an initial condition, all nodal temperatures are assumed to start with the same value as 
in the case when the Shuttle bay doors are just opened. In this analysis, the initial nodal 
temperature was taken to be 250 K. However, it is only the base that is assumed to have a 
large enough thermal mass so as to remain essentially constant at this temperature 
throughout the simulation. The output of this program is a time vector and a temperature 
matrix which provides a temperature for each node at each instant in time in the 
integrating process. For more details of the program implementation, reference may be 
made to the listing of the Matlab function which is provided in Appendix D. The user 
inputs the initial start time (to), the final stop time (tf) in seconds, and the uniform initial 


temperature in a 12 X 1 initial condition vector (TO) in this Matlab function. 
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the “‘ode45” integration procedure is carried out with the help of this function to integrate 


the coupled equations in time. 


Several representative cases were run to see the trends in nodal temperatures for 
different values of the solar incidence angle and the strip heater power input. The 


following cases were tried: 


SOLAR ~—__| HEATER POWER - 


ANGLE 
(DEGREES) 





Table 4.1. Test cases used for analysis (in the canister) 


It may be recalled that each heater has a 22 W capacity. All the heaters were either on or 
off during the simulation. Each provided an even distribution of heat to the lower deck of 
PANSAT. The other components of PANSAT are assumed to have no power dissipation 
in these cases. 


2. PANSAT in orbit 


The same methodology is the same as was used above in part 1, except that there 
are eleven nodes, and hence eleven nodal equations. The base (node 11) temperature is 


not fixed to the base temperature as in the previous case and are allowed to evolve 
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naturally. Table 4.2 presents the cases that were viewed for the in-orbit case. The 
orientations are shown in Figure 14 below and are assumed to be constant over the entire 


orbit. 


CASE# | SOLAR FLUX | ROTATION “HEAT ORIENTATION | 
RATE an ee 


Wo anti-earth 







4A ioee = = node | 
| 4B 1309 | node | 


SA _| 1399 i 
ee _).1309 ee eee Sica 


Table 4.2. Test cases used for analysis (in orbit) 
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Figure 14. Orientations of PANSAT 
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VI. RESULTS AND DISCUSSION 


The results and discussion are presented below for the two cases of PANSAT (A) 
the gas canister and (B) in orbit. The time integration is only for 90 minutes because the 
period of the orbit is approximately 90 minutes based on the scheduled altitude of 574 km 
as shown in Figure 3. The altitude of the mission could change. However, the space 
shuttle has maximum altitude constraints. Furthermore, PANSAT has no propulsion to 
change its altitude once in orbit. Hence; PANSAT should have a period of approximately 
90 minutes with an eclipse time of approximately 30 minutes throughout its lifetime. 

A typical low earth orbit has a period of approximately 90 minutes and an eclipse 
period of approximately 30 minutes. The period increases with an increase in altitude. 


Conversely, the eclipse time decreases as the altitude increases. 


A. PANSAT IN THE GAS CANISTER 


The results from the six cases tabulated above are depicted in the following pages 
in the form of temperature-time plots showing their approach to steady state conditions. 
The first graph in each set shows the variation of the temperature of the external nodes 
(1,6,7,8,11,12) with time. The second graph of each set is the temperature of the internal 


nodes (2,3,4,5,9,10) as it varies with time. 
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Figure 15. Nodal temperatures vs time. case 1A. 


Figure 15 shows the temperature variation for normal solar incidence (zero angle) 


with the heaters off. Note that nodes 11 and 12 (the GAS canister walls) experience the 
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lowest temperatures since they do not receive any sunlight in this orientation, and are 
exposed to cold black space. The highest temperature is felt at node | because it receives 
direct sunlight at normal incidence. The temperature of node 6 does not change 
appreciably. There are several reasons for this behavior: node 6 has several conduction 
paths to lose heat (via nodes 7,8,11,12,base). In particular the conduction between node 6 
and the base of the GAS canister is very large due to the large thermal mass of the base 
and its role as an effective heat sink. Thus there is a short circuit effect between node 6 
and the base which plays a dominant role in the steady state and transient analyses. 
Furthermore, the internal nodes (2,3,4,5,9,10) are bounded by the temperature of node 1 


and node 6. 
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Figure 16. Nodal temperatures vs time. case 1B 
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Figure 16 is the same case as in Figure 15 except that the heaters are turned on. Recall 
that the combined power of the heaters is 88 W. This power has a negligible heating 
effect on node 6 (and PANSAT) primarily due to the short-circuit effect noted earlier by 
which heat transferred to node 6 is more easily conducted away through the base, rather 


than being transferred into PANSAT. 
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Figure 17. Nodal temperatures vs time. case 2A 
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Figure 17 depicts the temperature variation with time for a solar incidence angle 
of 45 degrees with the heaters off. The net effect of this new configuration is that the heat 
input to node | (top panels) is reduced, while the heat input to node 11 (one half of the 
GAS canister) is increased. As a result the temperature of node 1 is less than that in the 
normal incidence case . The temperature of node 11 is higher than that of node 12 since it 
receives direct sunlight whereas node 12 does not. Also, the temperature of node 7 is 
slightly greater than node 8. In contrast, note that in Figures 12 and 13, nodes 11 and 12 
were at roughly the same temperature. Once again, the temperature of node 6 remains 


relatively constant over the time period it takes to reach steady state. 
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Figure 18. Nodal temperatures vs time. case 2B. 
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As discussed before with regard to Figure 17 the status of the heaters is fairly 
inconsequential to the problem, and the temperature-time history is not very different 


from that in Figure 17. The same reasons are applicable. 
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Figure 19. Nodal temperatures vs time. case 3A. 
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In this case the solar angle of incidence is 90 degrees, i.e. node I receives no 
direct sunlight, while node 11 receives all of the direct sunlight. Consequently the 
temperature of node 11 is higher, while that of node | is lower, than that in Figure 16 for 
the 45 degree incidence case. Here too the temperature of node 6 remains largely 


constant. 
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Figure 20. Nodal temperatures vs time. case 3B. 


The temperature variation in Figure 20 is similar to Figure 19. The strip heater 


power has no significant effect on the temperature behavior. The large conductance 
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between the base and node 6 dominates the energy exchange as discussed earlier. Again 
the temperature of node 1 and node 6 bound the temperatures all other nodes of 


PANSAT. 


B. PANSAT IN ORBIT 


The test cases for PANSAT on orbit are presented in the temperature-time plots 
below. The first graph shows the temperature external nodes (1,6,7,8,11). The second 
graph shows the temperature of the internal nodes (2,3,4,5,9,10). The numbers 
corresponding to the nodal temperature is shown in each figure. The period for 
PANSAT’s orbit is approximately 90 minutes with and eclipse time of 30 minutes. 

Also shown are temperature-time plots of the batteries over an extended period 
(approximately 30 orbits) which show the evolution of the temperature variation to its 


final periodic state. 


49 


temperatur 





265 


DEW hese = 22 } Ot eee 





temperature(K) 


2a fon ne enn =e 


—_ eae eeeee Freee eae @& 


250 


time(hrs) 


Figure 21. Nodal temperatures vs time. case 4A. 


Figure 21 depicts an orientation with node 7 facing the sun. Note that the temperature of 
node 1 and node 7 are approximately the same value. Nodes 1, 6, and 11 receive both 


ene 


solar and earth emitted radiation. The temperature of node 1 is eclipsed by the 
temperature of node 7 for approximately 30 minutes until the effect of the earth emitted 
radiation causes the temperature of node 1 to surpass the temperature of node 7. Nodes 6 


and 11 respond the slowest to temperature gradients in this arrangement. 
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Figure 22. Nodal tempeartures vs time. case 4B 


a2 


Figure 22 depicts the same general orientation as in Figure 21. PANSAT has a rotation 
rate of 4 degrees per second. The minimum solar constant value is used and there is no 
internal heat dissipation. The rotation of the spacecraft minimizes the development of 
temperature gradients. Hence, the external surfaces of PANSAT radiate heat to cold 
black space. Again, the time phase lag is evident. The temperature of Node 1 decreases 


most rapidly. The nodal temperatures are bounded by node | and node 6. 
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Figure 23. Nodal temperatures vs time. case 5A 
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Figure 23 depicts the orientation with node 1 facing the sun. Nodes 7 and 8 have 
approximately equal temperatures values. Since node | receives direct sunlight, it reacts 
more quickly to the temperature gradient. Each of the remaining nodal temperatures 
exhibit a time-phase lag. The greater the thermal distance between a given node and node 
1, the greater the time-phase lag. Thus, node 6 and node 11 have the most pronounced 


time-phase lag. The temperature of node 1 and 6 acts as a boundary for PANSAT. 
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Figure 24. Nodal temperatures vs time. case 5B 


Figure 24 depicts the PANSAT with a rotation rate of 4 degrees per second. The results 


are similar to Figure 22 and the same reasons apply to this case also. 
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Figure 25. Temperature of Battery A and Battery B vs time. Case 4A 
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Figure 25 shows the variation of the temperature of the batteries (A and B) over a period 
of 2 days. In both plots the temperatures reach steady state values in approximately 18 


hours and oscillate between 281 K and 286 K with an average temperature of 284 K. 
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Figure 26. Temperature of Battery A and Battery B vs time. Case 4A 
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Figure 26 shows the same temperature-time variation as in figure 25. The difference 
between Figure 25 and Figure 26 is that in this scenario there are no internal heat 
dissipation from the spacecraft electronic components. Hence, the steady state 


temperature values are approximately 15 K lower than in the previous case. 
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Figure 27. Temperature of battery A and battery A vs time. Case 5A 
Figure 27 depicts the temperature for the orientation in case 5A. Again, the temperatures 


stabililize in approximately 18 hours at a temperature of 278 K. 
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Figure 28. Temperature of battery A and battery A vs time. Case 5A 
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Figure 28 shows the same scenario as Figure 27 does except that there is no internal heat 
dissipation in this scenario. A slightly larger temperature difference is observed when 
there is no internal heat dissipation as compare to the case examined in Figures 25 and 


2G. 
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VII. PANSAT THERMAL TESTING 


A. THERMAL TESTING IN THE DESIGN PROCESS 

Thermal testing is used to verify the design process and increase the confidence 
in the overall design. System as well as subsystem level testing is conducted. Thermal 
design is iterative process requiring analysis to improve the thermal performance of the 
spacecraft. The analytical model is subsequently modified to accurately represent the 
thermal environment and thermal boundaries. 

MIL-STD-1540 B and MIL-STD-343, reference guides for system and component 
testing, details testing requirements for military spacecraft [Ref. 1]. All test plans will 
incorporate some elements of this guide. Some tests may not be appropriate for all 


spacecrafts. 


B. SYSTEM LEVEL THERMAL TESTING 

A test plan has been developed for PANSAT. It is scheduled to undergo a thermal 
balance test to validate the thermal analytical model. The thermal balance test will be 
conducted in conjunction with the thermal vacuum test. This test would normally 
identify major design errors in the thermal control system [Ref. 7]. Each of the electronic 
components will undergo thermal testing. The thermal vacuum test verify the proper 
interaction of components and subsystems. 

The number and type of tests conducted vary from spacecraft to spacecraft. 
Determining the right assortment of tests takes into account the mission’s objective, 


costs, facility availability and schedule. 
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VIII. CONCLUSIONS 


As stated in Chapter I, the emphasis of this analysis was to determine 
temperature-time history of PANSAT and the temperatures of other key components in 
PANSAT such as the two Ni-Cd battery compartments. The goal of this work is to 
ascertain that these temperatures are maintained within their design limits while 
PANSAT is in a stowed condition awaiting launch and while in orbit. Based on the 
requirements and assumptions presented in Chapters IJ and IV respectively, the following 


conclusions may be made: 


A. PANSAT IN THE CANISTER 


e In all six cases presented in Chapter V, the battery compartments temperature remain 


well within the required limit (-30C to +60C). 


e The worst scenario for the battery compartments is in cases 3A and 3B for which the 
solar angle of incidence is 90 degrees and the compartments receive more heat than 


for the case of normal incidence from directly overhead. 


e The base of the GAS canister is thermally large and serves effectively as a heat sink 
which dominates the thermal behavior of the system by providing a low resistance 


heat flow path from PANSAT to the rest of the Shuttle. 


e The battery compartments exchange more heat by conduction with the lower 
equipment plate (node 5) than by radiation with the digital control system (node 4) 


and the mid-deck panels (nodes 7 and 8). 


e The kapton strip heaters have a negligible impact on the thermal behavior of 
PANSAT since the heat that they radiate to node 6 is more easily conducted away 


through the base, rather than transferred into PANSAT. 


e The difference between the steady state values and the initial value of the 


temperatures has a weak dependence on the choice of the initial temperature. This is 


ow 


because the absolute value of the temperature determines only the radiation 
conductances which are relatively small. Hence, a different initial temperature 
effectively results in a corresponding up-down shift in the temperature-time plots 


discussed earlier. 


e Considering the low temperatures experienced by PANSAT while it is stowed in the 
GAS canister, it’d be desirable to expose PANSAT to a “sun-soak”’ prior to launch to 


ensure that the temperatures are closer to the operating values. 


e Such a model and the attached program can also be used to simulate cyclic 
heating/cooling conditions such as those experienced by PANSAT while in orbit. 
With some modifications the program can also be made to account for the variation of 


the solar incidence angle (and hence the solar heat input) with time. 


B. PANSAT IN ORBIT 


The following conclusions can be made for the in-orbit phase of the 

mission. 

e The temperature of node 1 is the highest and the temperature of node 6 
is the lowest in both orientations. All other nodal temperatures are 
within these temperature limits. 

e The time phase lag of the temperature response of the nodes verify the 
accuracy of the model. Nodes (6,11), the furthermost distance from 
node | in exhibit the largest time phase lag and amplitude dampening. 
Node 11 reacts more slowly than node 6 due to its much larger mass. 

e The battery compartment temperatures remain within the specified 
temperature limits (-5 C to 35 C) for all orientations with no rotation 


Tae. 
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This procedure is not dependent on the initial temperature of the 
nodes. The initial nodal temperatures affect the rise-time to steady state 
and not the final steady state temperature values. 

Further testing needs to be completed to determine the temperature 
profile for cases with different rotation rates. 

The temperature profiles of cases 4B and 5B, PANSAT with a rotation 
rate, exhibit the lowest absolute temperature margins of all the cases. 
However, the temperature differences is consistent with all of the other 
cases. 

PANSAT should receive a sun-soak prior to insertion into its orbit to 
ensure that the initial nodal temperatures are stabilized near normal 
operating values. This would reduce the time required to bring 


PANSAT to full operations. 
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APPENDIX A. NODAL DESCRIPTION 


From To Conductance [W/K] Nodal Description 














For conduction from the top and upper deck panels, and the 


triangular antenna subassembly panel to the mid-deck panels 


For conduction from the top and upper deck panels, and the 


triangular antenna subassembly panels to the mid-deck panels 


— 





oe - a, —<—— ee 


For conduction from the mid-deck panels to the bottom deck panel 
and lower deck solar panels and the triangular antenna 


| subassembly panels 


For conduction from the mid-deck panels to the bottom deck panel 
and the lower deck panels and the 4 triangular antenna 


subassembly panels 


For radiation from top and upper deck panels and triangular 





antenna subassembly panels to the mass storage containers A and 


B and multiplexer 


2 3 | C,,=126 For conduction from mass storage and multiplexer to the upper 


equipment plate. 





Bees @| C,.=0-26 For radiation from the upper equipment plate and the electrical 


power system to the digital control system 





4 5/7 C,,=108 | For conduction from the digital control system to lower 


equipment plate and the radio frequency communications 


Weil 


ee For radiation from battery A to digital control system 


8 










|For conduction from the lower equipment plate and radio | 





= sD 


frequency communications to lower deck panels and the bottom 


panel 


For conduction from the mid-deck panels to the lower equipment 


plate and radio frequency communications 


For transfer by conduction from the mid-deck panels to the lower 


equipment plate and radio frequency communications 


For radiation from the mid-deck panels to Battery A 


Sia For radiation from battery B to The digital control system 


Co, 
ST pe 
Se ilic 






=92 
=29.8 
= 29 










For radiation from the mid-deck panels to the battery B 


For conduction from the mid-deck panels to upper equipment plate | 


and electrical power system(EPS) 


For conduction from the 4 mid-deck panels to Upper equipment 





plate and the electrical power system 


For conduction from the bottom panel through the Launch vehicle 


interface to the base of the GAS canister 


For conductance from battery B to the lower equipment plate and 


the radio frequency communications 


For conductance from battery A to the lower equipment plate and 





the radio frequency communications 


az 


APPENDIX B. NODAL COMPOSITION 
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top panel, 4 upper deck panels ,4 triangular 
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upper equipment plate, electrical power system 
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bottom panel, 4 lower deck panels, 4 triangular | 3.08 
panels 


4 mid-deck panels 


8 _| 4 mid-deck panels 163 «| 0.07 __ 
9 | battery A | 1.97 0.02 


launch vehicle interface 
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APPENDIX C. PROGRAM CODE 


% The following code describes a Matlab function used to carry out a transient thermal 
analysis of the temperatures of a nodal thermal model of PANSAT. This program is used 
for the cases 1A-3A. (PANSAT in the GAS canister) 


function Tdot=thm0Ol(t, T) 
% The inputs ‘t’ and’ T’ represent the time (to, tf) and the initial nodal temperature 


%oV%o%%% ail temperatures are indegreesK %S%UMhMHHLLALAKHLKL% 
% The temperature of the nth node is represented as T(n), n=1..12 


% "sig" below is the Stefan-Boltzmann constant [W/m‘2-K‘4] 
sig= 5.67e-8; 


% emissivity and absorptivity of aluminum (dimensionless) 
% assume mass of solar cells much less than underlying aluminum panel 


em=0.83; % emissivity of solar cells 
embc=0.88; % emissivity of beta cloth 
emal=0.80; % emissivity of aluminum 
abs=0.78; % absorptivity of solar cells 
absbc=0.30; % absorptivity of beta cloth 
absal=0.30; % absorptivity of aluminum 


% The F's below are the view factors between the different nodes (dimensionless) 


F12=0.16; %assume F21=1 
F34=0.63; Yassume F43=1 
F79=0.23; %assume F97=1 
F94=0.63; %assume F49=]1 
F810=0.23; Z%assume F108=1 
F104=0.63; %assume F410=1 
F811=1; 

Ey 12=1. 

Fh6=0.03; 

F6c=0.38; 


% masses of nodes [kg] 


i> 


ml=3.08;: % mass of node | 


mic=2. 18: 
Mio=2.6 ||; 
m4=4.14; 
mo=3.67 
m6=—s 0s, 
mas .05; 
m&=1.63; 
m9=1-9 7: 


miO=197/-: 
fil). 29- 
mi 2=1.293 


%%%% specific heat capacity of aluminum [J/kg-K] %%%D% 


Cp=903; 


% The A’s below are the nodal heat transfer areas [m2] 


Ad=O 15: 
A6=0.15; 
A6b=0.04; 
A6oc=0.11; 
Ac6=0.60; 
Al=0.03; 
A3=0.06; 
A4=0.04; 
A7=0.07; 
A8=0.07; 
A9=0.02; 
A10=0.02; 
Al1=0.14; 
A12=0.14; 


% area of node 6 that sees the can 
% area of can that sees node 6 


% Radiation conductances in the thermal network [W/K] %2%%%%% 


C12= (sig*(T(1).42 + T(2).42)*(TU)+T(2))/((1-emal)/(emal* A1) +(1/(A1*F12)) +... 
(1-emal)/(emal*A2))); 


C34= (sig*(T(3).42+T(4).42)*(T(3)+T(4))/((1-emal)/(emal* A3)+(1/A3*F34)+(1-emal)/... 


(emal*A4))); 


C94=(sig*(T(9).42+T(4).42)*(T(9)+T(4))/((1-emal)/(emal *A9) +(1/A9*F94)+(1-emal)/ 


(emal * A4))); 


1G 


C104=(sig*(T(10).42+T(4).%2)*(T(10)+T(4))/((1-emal)/(emal* A 10)+(1/A 10*F104)+(1- 
emal)/(emal * A4))); 


C79= (sig*(T(7).42+T(9).42)*(T(7)+T(9))/(C1 -emal)/(emal* A7)+( 1/A7*F79)... 
+(1-emal)/(emal*A9))) ; 


C8 10= (sig*(T(8).42+T(10).%2)*(T(8)+T(10))/((1-emal)/(emal* A8)+(1/A8*F 810)... 
+(1-emal)/(emal*A10))); 


% Conduction conductances in the thermal network [W/K] 


©23=1206; % node 2 to 3 
@is—2.10; % node 1 to 8 
el —=2.10; % node | to 7 
C86=2.10; % node 8 to 6 
@€76—2.10;: % node 7 to 6 
7 5—3.50; % node 7 to 5 
Co5=5.0; % node 8 to 3 
9 3=3.50; % node 7 to 3 
@55=3.50; % node 8 to 5 
C56=0.98; % node 5 to 6 
C45=108; % node 10 to 8 
Cé6b=92; % node 6 to base 
C95=29.8; % node 9 to 5 
C105=29.8; %node 1l0to5 


% Temperature declarations {K] 2%% 


Tinit=250; 9% initial temperature of PANSAT prior to opening the Shuttle Hatch 
T_space =0; % temperature of space 


Tc6 =.5.*(T(11)+T(12)); %temperature of the part of the canister that sees node6 
assumed to be an average of T911) and T(12). 


Tb=Tinit; % Temperature of the base of the shuttle compartment. 
% The mass of the base is much larger lvi 
ToTo TPT %%%% Heat Input(W] 26MMOLDLUMHHHLM% 


qsol=1353; % average Solar constant [W/m‘2] 

theta=pi/2; % angle of solar incidence (90 degrees in this example) 
Qh=88; % input power from all 4 strip heaters to node 6 
Qh6in= Qh*abs* Fh6; % heater input directed to node 6 


Q1in=qsol*cos(theta)* A 1 *abs; 
Q1 lin=qsol*sin(theta)*A 1 1*absbc; 


To 


Q12in=0; 


Too Vo To To FVo%O%%%o%o Heat Loss [(W] %2%%2LLDLD%% 


Q6b=C6b.*(T(6)-Tb); 

Qlsp=(em.*sig.*(T(1).%4-T_space.*4)*A 1); 
Qb6c=(sig.*(T(6).*4-Tc6.%4))./((1-em)/(em*A6c) + 1/(A6c*F6c)+(1-emal)/(emal* Ac6)); 
Q7c=(sig.*(T(7).%4-T(12).%4)).((1-em)/(em*A7) + 1//A7*F712)+(1-emal)/(emal*A 12)); 
O8&c=(sig.*(T(8).%4-TC1 1).%4))./((1-em)/(em*A8) + 1/(A8*F811)+(1-emal)/(emal*A11)); 
Q11lsp=embc.* sig.* (T(11).%4-T_space.%4); 

Q12sp=embc.* sig.* (T(12).%4-T_space.“4); 


%% the differential equations from nodal energy balances %%%%% 


Tdotl= -(((C17+C184+C12).*T(1)-C12.*T(2)-C17.*T(7)-C18.*T(8)- 
lin+Q1sp)/(m1*Cp)); 

Tdot2= -((-C12.*T(1)+(C12+C23).*T(2)-C23.*T(3))/(m2*Cp)); 
Tdot3=-((-C23.*T(2)+(C73+C23+C83+C34).*T(3)-C34.*T(4)-C73.*T(7)- 
...C83.*T(8))/(m3*Cp)); 
Tdot4=-((-C34.*T(3)+(C34+C45+C104+C94).*T(4)-C45.*T(5)-C94.*T(9)- 
...C104.*T(10))/(m4*Cp)); 
Tdot5=-((-C45.*T(4)+(C75+C45+C85+C56).*T(5)-C56.*T(6)-C75.*T(7)-C85.*T(8)- 
C95.* ...(T(9)-T(S))-C105.*(T(10)-T(5)))/(m5*Cp)); 
Tdot6=-((-C76.*T(7)-C86.*T(8)-C56.*T(5)+(C56+C86+C76).*T(6)+Q6b+Q6c- 
..-Qh6in)/(m6*Cp)); 

Tdot7=-((-C17.*T(1)-C73.*T(3)-C75.*T(5)- 
C76.*T(6)+(C76+C17+C734+C79+C75).*T(7)-C79.* T(9)+Q7c)/(m7*Cp)); 
Tdot8=-((-C18.*T(1)-C83.*T(3)-C85.*T(5)-C86.*T(6)- 
C810.*T(10)+(C85+C18+C83+C86+C8 10) .*T(8)+Q8c)/(m8*Cp)); 
Tdot9=-((-C94.*T(4)-C79.*T(7)+(C79+C94).*T(9)+C95.*T(9)-C95.*T(5))/(m9*Cp)); 
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Tdotl10=-((-C104.*T(4)-C810.*T(8)+ (C8 10+C104).*T(10)+C105.*T(10)- 
meena. 1(5))/Gnl0*Cp)): 

Tdot! 1=-((-Q1 1in+Q1 1sp-Q8c)/(m11*Cp)); 

Tdot 12=-((-Q12in+Q12sp-Q7c)/(m12*Cp)); 


%%I%%% now assemble all the Tdot’s into one vector Tdot 2%™%%MHM% 


Tdot=[Tdot! Tdot2 Tdot3 Tdot4 Tdot5S Tdot6 Tdot7 Tdot8 Tdot9 TdotlO Tdot11 
Tdot12]; 


% The time integration procedure is now carried out with Tdot as input to the “ode45” 


routine in Matlab along with the start and finish times. 


ag 


% This program is used for the cases 4A and 4B. 
function Tdot=thm21vs(t, T) 
%%all temperatures are in degrees K %%%%%V%%o%o Vo Vo Vo Vo Vo Vo 


% "sig" is the Stefan-Boltzmann constant (mks units: W/m42-K‘%4) 
sig= 5.67e-8; 


% emissivity and absorptivity of aluminum (dimensionless) 
% assume mass of solar panel much less than underlying panel 


eMm=0 sa: 
emal=0.80; 
abs=0.78; 


_ absal=0.30; 


% F’s are the view factors (dimensionless) 
F12=0.16; %assume F21=1 

F34=0.63; Yassume F43=1 

F79=0.23; Yassume F97=1 

F94=0.63; Yassume F49=1 

F810=0.23; Zassume F108=1 
F104=0.63; Zassume F410=1 


% masses of nodes (Kg) 


m1=3.08; % mass of node | 
m2=2.18; 
m3=261; 
m4=4.14; 
Mo=3°07 
m6=3.08; 
m7=1.63; 
m8&=1.63: 
m9=1.97; 
m10=1.97; 
m11=5.68; 


%o%o Yospecific heat capacityJ/kg-K)%%%%%%o% ToT Mo To Vo Yo To Yo Yo Yo Yo Yo To To To Yo To 
Cp=903; 


% A’s are the nodal areas (in m‘2) 
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exit =0.15; 
A6=0.15; 
A6b=0.04; 
A2=0.03; 
A3=0.06; 
A4=0.04; 
A7=0.07; 
A8=0.07; 
A9=0.02; 
A10=0.02; 
per =0.17; 


%%% radiation conductances in the thermal network (W/K)%%%%%% 
% The Temperatures are the Temperatures at a given time "t" % 


C12= (sig*(T(1).42 + T(2).42)*(T(1)+T(2))/((1-emal)/(emal*A 1) +(1/(A1*F12)) 
+(1-emal)/(emal* A2))); 


C34= (sig*(T(3).42+T(4).%2)*(1(3)+T(4))/(C1-emal)/(emal* A3)+(1/A3*F34)+(1-emal)/... 
(emal*A4))) 


C94= (sig*(T(Q).424+T(4).42)*(T9)+T(4))/((1-emal)/(emal * A9) +(1/A9*F94)+(1- 
emal)/(emal *A4))); 


C104=(sig*(T(10).42+T(4).%2)*(T(10)+T(4))/((1-emal)(emal* A10)+(1/A10*F104)+(1- 
emal)/(emal * A4))); 


C79= (sig*(T(7).42+T(9).42)*(T(7)+T(9))/((1-emal)/(emal* A7)+(1/A7*F79)... 
+(1-emal)/(emal* A9))) ; 


C810=(sig*(T(8).42+T(10).2)*(T(8)+T(10))/((1-emal)/(emal* A8)+(1/A8*F8 10)... 
+(1-emal)/(emal*A10))); 


% Conduction Conductances in the thermal network (W/K) %%%%% 


C23=126; % node 2 to 3 
C18=2.10; % node 1 to 8 
C17=2.10; % node | to 7 
C86=2.10; % node 8 to 6 
C76=2.10; % node 7 to 6 
C75=3.50; % node 7 to 5 
C83=3.50; % node 8 to 3 
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C73=3.50; % node 7 to 3 
C85=3.50:; % node & to 5 
C56=0.98; % node 5 to 6 
C45=108; % node 10 to 8 
C6b=92; % node 6 to base 
C95=29.8; % node 9 to5 
C105=29.8; % node 10 to 5 


% %%oTemperature intializations (in K) %%%%%%%V%%o%oV%o Fo Vo Mo 
Tinit=250; %initial temperature of the Pansat prior to opening the Shuttle Hatch 


T_space = 0; % Temperature of space 


JoTo Too %o%o Mo Heat Input (in W) %%%%%%%o% Vo To To Yo To Mo To To Yo Yo Yo Vo Yo Yo Vo To Yo Yo 
omega=0; 


if t< 3600 % period is approximately 1.5 hrs with 
% an eclipse time of approximately 
Jo 1800sec 
iis@l—b 370. 
q_albedo=77.1; 


else 


qsol=0; 
q_albedo=0; 


end 


theta=0; % solar angle of incidence 
q_earth=32.08; % earth emitted radiation 


Qlin=qsol.* cos(omega.*t)*cos(theta)*A 1 *abs; 
Q6in=(q_earth +q_albedo)*A6*abs; 

© in=0: 

Qsin=0; 

Q1lin=(q_earth +q_albedo)*(A6b)*absal; 
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ToT To Vo To Yo Yo Yo Fo YoVoequipment power dissipations-W %% % % % Yo Yo Yo To Go Yo Yo Vo 


q_eps=0; % electrical power system 

q_bata=0; %battery a 

q_batb=0; % battery b 

q_rfcomm=0;% Communication system (receiver and transmitter) 
q_dcs=0; % digital control system 


%o%o %o%o Vo % Heat Loss (in W) %2%%%%%% Go % No Vo Go Mo Yo Yo Yo Yo Yo To Yo Yo Yo Ge Vo Yo Yo 


Q6b=C6b.*(T(6)-T(11)); 
Q1sp=(em.*sig.*(T(1).*4-T_space.*4)*A1); 
Qb6sp=(em.*sig.*(T(6).%4-T_space.*4)*A6); 
Q7sp=(em.*sig.*(T(7).%4-T_space.*4)*A7); 
OQ8sp=(em.*sig.*(T(8).4-T_space.*4)*A8); 

Q1 1sp=(emal.*sig.*(T(11).44-T_space.*4)*A6b); 


%%o%o%% the differential equations %%%%%% 7% %%o Yoo Yo Yo Yo To Yo Yo Yo Yo Vo 


Tdotl= -(((C17+C18+C12).*T(1)-C12.*T(2)-C17.*T(7)-C18.*T(8)- 
Qlin+Q1sp)/m1*Cp)); 


Tdot2= -((-C12.*T(1)+(C12+C23).*T(2)-C23.*T(3))/(m2*Cp)); 


Tdot3=-((-C23. Re eae *T(3)-C34.*T(4)-C73.*T(7)-C83.*T(8)- 
q_eps)/(m3*Cp)); 


Tdot4= -((-C34.*T(3)+(C344+C454+C104+C94).*T(4)-C45.*T(5)-C94.*T(9)- 
C104.*T(10)-q_dcs)/(m4*Cp)); 


Tdot5= -((-C45.*T(4)+(C75+C45+C85+C56).*T(5)-C56.*T(6)-C75.*T(7)-C85.*T(8)- 
C95.* (T(9)-T(S))-C105.*(T(10)-T(5))-q_rfcomm)/(m5*Cp)); 


Tdot6= -((-C76.*T(7)-C86.*T(8)- C56.*T(5)+(C56+C86+C76).*T(6)+Q6b+Q6sp- 
Q6in)/(m6*Cp)); 


Tdot7=-((-C17.*T(1)-C73.*T(3)-C75.*T(5)-C76.*T(6)+... 
(C76+C 17+C73+C79+C75).*T(7)-C79.*T(9)-Q7in+Q7sp)/(m7*Cp)); 


Tdot8=-((-C18.*T(1)-C83.*T(3)-C85.*T(5)-C86.*T(6)- 
C810.*T(10)+(C85+C18+C83+C86+C8 10).*T(8)-Q8in+Q8sp)/(m8*Cp)): 
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Tdot9=-((-C94.*T(4)-C79.*T(7)4+(C794+C94).*T(9)+C95.*T(9)-C95.*T(5)- 
q_bata)/(m9*Cp)); 


Tdot10=-((-C 104.*T(4)-C8 10. oe (C8 10+C104).*T(10)+C105.*T(10)-C105.*T(5)- 
q_batb)/(m10*Cp)); 


Tdot 1 1=-((-Q6b+Q1 1sp-Q1 1in)(m11*Cp)): 


%%o%% now assemble all the "Tdot"s %%%%%%%% %o%o Go Vo Vo Vo To Fo Vo To Yo Vo 


Tdot = [Tdot! Tdot2 Tdot3 Tdot4 Tdot5 Tdot6 Tdot7 Tdot8 Tdot9 Tdot10 Tdot1 1]; 


function Tdot=thm11vs(t, T) 
% This program is used for cases 5A and 5B 


To To To To Too Yo Yo Yo Yo %o%o all temperatures are in degrees K %%%%%%%%%oFo% Vo To Vo 


% "sig" is the Stefan-Boltzmann constant (mks units: W/m42-K‘4) 
sig= 5.67e-8; 


J emissivity and absorptivity of aluminum (dimensionless) 

% assume mass of solar panel much less than underlying panel? 
em=0.83; 

emal=0.80; 

abs=0.78; 

absal=0.30; 


% F’s are the view factors (dimensionless) 
F12=0.16; Y%assume F21=1 

F34=0.63; %assume F43=1 

F/79=0.23; Yassume F97=1 

F94=0.63; %assume F49=1 

F8 10=0.23; Yoassume F108=1 

F104=0.63; %assume F410=1 
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% masses of nodes (Kg) 


m1=3.08; % mass of node | 
m2=2.18: 
me=2.61; 
m4=4.14; 
Mmo—=o.0/ ; 
m6=3.08; 
m= 1ho3: 
Mmio=!.03; 
m9=1.97; 
m10=1.97; 
m11=5.68; 


J%o To To Vo Yo Yo Yospecific heat capacity(J/kg-K)%%%% % To % To To Yo Yo M0 Yo M0 Yo Vo To Yo Yo Yo 
Cp=903; 
% A’s are the nodal areas (in m“2) 


mal=0.15; 
Alp=0.06; - 
A6=0.15; 
A6b=0.04; 
A6p=0.06; 
A2=0.03; 
A3=0.06; 
A4=0.04; 
A7=0.07; 
A8=0.07; 
A9=0.02; 
A10=0.02; 
A11=0.04; 
Al11p=0.01; 


Jo To %o Yo Yo Yo Yo Radiation Conductances in the thermal network (W/K) %%%%%%%% % 
% The Temperatures are the Temperatures at a given time "t” 

Jo 

C12= (sig*(T(1).42 + T(2).42)*(T(1)+T(2))/((1-emal)/(emal*A1) +(1/(A1*F12)) +... 
(1-emal)/(emal* A2))); 
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C34= (sig*(T(3).*2+T(4).%2)*(T(3)+1(4))/((1-emal)/(emal* A3)+(1/A3*F34)+(1-emal)/... 
(emal* A4))); 


C94= (sig*(T(9).424+T(4).%2)*(T(9)4+1(4))/((1-emal)/(emal *A9) +(1/A9*F94)+(1- 
emal)/...(emal *A4))); 


C104=(sig*(T(10).%2+T(4).%2)*(T(10)+T(4))/((1-emal)/(emal* A10)+(1/A10*F104)+(1- 
emal)/...cemal * A4))); 


C79= (sig*(T(7).42+T(9).42)*(T(7)4+T(9))/((1-emal)/(emal*A7)+(1/A7*F79)... 
+(1-emal)/(emal*A9))) ; 


C8 10= (sig*(T(8).*2+T(10).%2)*(1(8)+T(10))/((1-emal)/(emal* A8)+(1/A8*F8 10)... 
+(1-emal)/(emal*A10))); 


Jo%o%o%o% Conduction Conductances in the thermal network (W/K) 4%%%%MH2RM% 


C23=126; % node 2 to 3 
C18=2.10; % node 1 to 8 
C17=2.10; % node 1 to 7 
C86=2.10; % node 8 to 6 
C76=2.10; % node 7 to 6 
C75=3.50; % node 7 to 5 
C83=3.50; % node 8 to 3 
C73=3.50; % node 7 to 3 
C85=3.50; % node 8 to 5 
C56=0.98; % node 5 to 6 
C45=108; % node 10 to 8 
C6b=92; % node 6 to base 
C95=29.8; % node 9 to5 
C105=29.8; % node 10 to 5 


JoTo To To Vo Vo Vo VoVoVo%o Temperature intializations (in K) %%%%%%M%%%%VoVVo Vo 
Tinit=250; %initial temperature of the Pansat prior to departing the shuttle 

T_space = 0; % Temperature of space 

ToT Too % To%o FOV %o%o equipment power dissipations%%%h% %o% h%V% %o% To To Go Yo Yo 


q_eps=5; % electrical power system 
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q_bata=0; Zobattery a 

q_batb=5; % battery b 

q_rfcomm=8;% Communication system (receiver and transmitter) 
q_dcs=3; % digital control system 


Jo To Vo To Yo To Fo Yo Yo Heat Input (in W) %%%%%o%o To To Vo To To Vo Yo To To Yo To To Yo Yo To Vo Yo Vo 
omega=4*pi/180; % spin rate for pansat 
if t<3600 


qsol=1399; 
q_albedo=77.1; 


else 
qsol=0; 
q_albedo=0; 


end 
theta=90; 


q_earth=32.08; 


Q1in=qsol*cos(omega.*t)*sin(theta)*A 1p*abs+(q_earth+q_albedo)*A1p*abs; 
Q7in=qsol*cos(omega.*t)*sin(theta)*A7*abs; 

Q8in=(q_earth+q_albedo)*A8*abs; 

Q6in=qsol*cos(omega.*t)*sin(theta)* A6p*abs+(q_earth+q_albedo)* A6p*abs; 

Q1 lin=qsol*cos(omega.*t)*sin(theta)*absal* A 11p+(q_earth+q_albedo)*(A11p)*absal; 


JoTo To To To YoYo YoH eat Loss (in W)%%%% % % Vo To %o To To Vo To Vo To Yo To To Yo Yo Yo To Yo Vo Vo Yo 


Q6b=C6b.*(T(6)-T(11)); 
Qisp=(em.*sig.*(T(1).44-T_space.%4)*A 1); 
Qésp=(em.*sig.*(T(6).%4-T_space.*4)*A6); 
Q7sp=(em.*sig.*(T(7).44-T_space.*4)*A7); 
Q8sp=(em.*sig.*(T(8).%4-T_space.*4)*A8); 

Q1 lsp=(emal.*sig. *(T(1 1).44-T_space.*4)*A6b); 


ToToToPo PO %N%%% the differential equations %%%%%o% % Too Yo Vo Yo Yo Yo To Yo Yo Yo To Vo Vo 
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Tdotl= -(((C174+C18+C12).*TC1)-C12.*T(2)-C17.*T(7)-C18.* T(8)- 
Q1lin+QIsp)/(m1*Cp)); 


Mdot2==((-C 2.7 T(1)+(Ci2+C 23).* 1(2)-C23.* T(3))/Gm2*Cp)), 


Tdot3=-((-C23.*T(2)+(C73+C23+C83+C34).*T(3)-C34.*T(4)-C73.*T(7)-C83.*T(8)- 
Gueas)/(im3~Cp)), 


Tdot4= -((-C34.*T(3)+(C34+C45+C104+C94).*T(4)-C45.*T(5)-C94.*T(9)- 
C104.*T(10)-q_dces)/(m4*Cp)); 


TdotS= -((-C45.*T(4)+(C75+C45+C85+C56).*T(5)-C56.*T(6)-C75.*T(7)-C85.*T(8)- 
C95.* (T(9)-T(5))-C105.*(T(10)-T(5))-q_rfcomm)/(m5*Cp)); 

Tdot6= -((-C76.*T(7)-C86.*T(8)-C56.*T(5)+(C56+C86+C76).*T(6)+Q6b+Q6sp- 
Q6in)/(m6*Cp)); 


Tdot7=-((-C17.*T(1)-C73.*T(3)-C75.*T(5)-C76.*T(6)+... 
(C76+C 174+C73+C79+C75).*T(7)-C79.* T(9)-Q7in+Q7sp)/(m7*Cp)); 


Tdot8=-((-C18.*T(1)-C83.*T(3)-C85.*T(S)-C86.*T(6)-C8 10.*T(10)+... 
(C85+C18+C83+C86+C8 10).*T(8)-Q8in+Q8sp)/(m8*Cp)); 


Tdot9=-((-C94.*T(4)-C79.*T(7)+(C794+C94).*T(9)+C95.*T(9)-C95 .*T(5)- 
q_bata)/(m9*Cp)); 


Tdot10=-((-C104.*T(4)-C810.*T(8)+ (C8 104+C104).*T(10)+C105.*T(10)-C105.*T(5)- 
q_batb)/(m10*Cp)); 


Tdot1 1=-((-Q6b+Q 1 1sp-Q1 lin)/(m11*Cp)); 


To To %%o%o%%% now assemble all the “Tdot"s M%%%%h%V%o%% %o%o Vo % Yo Vo Mo Fo Yo Yo 


Tdot=[Tdotl Tdot2 Tdot3 Tdot4 TdotS Tdot6 Tdot7 Tdot8 Tdot9 Tdot10 Tdot11]; 
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